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I. INTRODUCTION 

Supersymmetric extensions of the Standard Model 
(SM) give natural framework for solving the hierarchy 
problem. They lead also to unification of interactions 
at energies of the order of mauT 10"^^ GeV. Super- 
symmetry (SUSY) itself is required by string theory, the 
best candidate by now which resolves most of the SM 
problems. 

The Minimal Supersymmetric Standard Model 
(MSSM) is the minimal in the interaction and particle 
content consistent extension of the SM. It provides 
rich phenomenology by the introduction of twice as 
many particles as we know from the SM, the so-called 
superpartners of usual particles. It possesses also an 
accidental symmetry called i?-parity, defined by means 
of the baryon {B), lepton (L), and spin {S) numbers 
as i? = (_i)3B-i-L-i-2S^ Conservation of i?-parity implies 
that SUSY particles cannot decay into non-SUSY ones. 
As a consequence the lightest SUSY particle must be 
stable and may be considered as a good candidate for 
dark matter. Theoretically, however, nothing motivates 
i?-parity conservation and many models which allow for 
breaking of this symmetry have been considered. Among 
them one can distinguish three groups. In the first one 
i?-parity violation (RpV) is introduced as a spontaneous 
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process triggered by a non-zero vacuum expectation 
value of some scalar field Other possibilities include 
the introduction of additional bi- |2| or trilinear 3] RpV 
terms in the superpotential. 

The most pressing motivation for looking on physics 
beyond the SM comes from the discovery of neutrino os- 
cillations 0, S IE ■ The experimental evidence for this 
process clearly indicates that neutrinos do have non-zero 
masses. What is more, the flavour and mass eigenstates 
are not the same, which leads to mixing between them. 
In experiments focused on neutrino oscillations one can 
measure the mixing angles as well as the differences of 
masses squared. For absolute values of masses one has 
to study the large scale structures of the Universe, the 
endpoint of the electron spectrum in beta decay of Tri- 
tium , or search for signals of the hypothetical neutri- 
noless double beta decay (OvfiP). The latter is also the 
only known process which distinguishes between Majo- 
rana and Dirac-like neutrinos. Its observation will prove 
the Majorana nature of these particles. 

The importance of theoretical studies of various prop- 
erties of neutrinos is obvious. In the present paper we 
focus on the transition magnetic moment, which is gener- 
ated by interactions of quark-squark and lepton-slepton 
self-energy loops with an external photon. This work is 
a natural continuation and extension of calculations pre- 
sented in Refs. 9, 1_0, 11, IJJ. 

In the following section we define the model, describe in 
detail our method of obtaining the supersymmetric par- 
ticle spectrum, and the way in which unification (GUT 
constraints) is imposed. In Sec. Ill we construct the 
neutrino mass matrices in loop mechanism, as well as cal- 
culate the transition magnetic moments. In our present 
approach the squark and quark mixings are taken into ac- 
count. We present our results and state the conclusions 
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in the last part of the paper. 



II. GUT CONSTRAINTS AND PARTICLE 
SPECTRUM 

Let us fix our framework to be the Minimal Super- 
symmetric Standard Model with supersymmetry break- 
ing mediated by the gravity force (SUGRA MSSM). We 
follow the conventions and notation from Ref. We 
allow for i?-parity violation by considering explicitly bi- 
linear and trilinear RpV terms in the superpotential. We 
focus, however, on the trilinear part, as the effects we 
want to study come from 1-loop processes, while the bi- 
linear terms induce neutrino masses at tree level. 

The i?-parity conserving part of the superpotential of 
MSSM has the form 



W 



MSSM 



while its RpV part reads 



(1) 



+ -exy.K.kUtDyDl + e^bn^L^Hl (2) 



The Y's are 3x3 Yukawa matrices. L and Q are the 
SU{2) left-handed doublets while E, U and D denote the 
right-handed lepton, up-quark and down-quark SU{2) 
singlets, respectively. Hi and H2 mean two Higgs dou- 
blets. We have introduced color indices x,y,z — 1,2,3, 
generation indices i,j,k — 1,2,3 and the SU(2) spinor 
indices a, 6, c = 1, 2. 

The introduction of RpV terms implies the existence of 
lepton or barion number violating processes, like the un- 
observed proton decay. Fortunatelly one may keep only 
one type of terms and it is not necessary to have both at 
one time. In order to get rid of too rapid proton decay 
and to allow for lepton number violating processes, like 
the neutrinoless double beta decay, it is customary to set 
A" = 0. 

We supply the model with scalar mass term 

= m^^hlhi+ml^h^+q^mlq + l^mil 
+ umfju^ + dm'jyd'^ + em|e^ (3) 

soft gauginos mass term {a — 1, 8 for gluinos) 
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^gaug. ^ ^ I^M^B^ B + NhWi W"- + Ahgjg^ + h.c. 

(4) 

as well as the supergravity mechanism of supersymmetry 
breaking, by introducing the Lagrangian 

£^°^* = eab[{AEhlth\ej + {Ao^qrhld,, 

+ {Auhqrh'^ujx + Biihlhl + B2ed1hl], (5) 



where lowercase letters stand for scalar components of 
respective chiral superfields, and 3x3 matrices A as well 
as B/i and B2 are the soft breaking coupling constants. 

The procedure of finding a GUT constrained low en- 
ergy spectrum of the model consists of a few steps. 

We take into account mass thresholds where SUSY 
particles start to contribute and use the SM 1-loop 
renormalization group equations (RGE) below ap- 
propriate threshold and 1-loop MSSM RGE [3 above 
it. The 2-loop corrections as well as corrections coming 
from the presence of RpV couplings are small and do not 
change the results by more than few percent ^3 • 

Initially all the thresholds are set to 1 TeV and are dy- 
namically modified during the running of mass parame- 
ters. The values of Yukawa couplings Y at mz are given 
by lepton and quark mass matrices M |l8j | 



Mu 
Md 
Me 



(6) 



where Vd = {Hi) and Vu — {H2) are the neutral Higgs 
vacuum expectation values. Their ratio defines the angle 
[3 through the relation Vu/vd — tan/3. S matrices per- 
form diagonalization so that one obtains eigenstates in 
the mass representation. 

After evolving the dimensionless couplings from the 
electroweak scale mz up to itigut = 2 x 10^^ GeV, we 
unify the masses of gauginos, sfermions and squarks to 
be equal to a common mass parameter toq. We set also 
the masses of all fermions to a common value mi/2- The 
trilinear soft couplings at mcuT are set according to the 
formula 

A, = AoYi, 

with Aq being another input parameter. We construct 
the squark, slepton, chargino and neutralino mass matri- 
ces and in the next step evolve all the quantities down to 
mz- During that running the tree-level Higgs potential 
is minimized, ie. the following set of equations is solved: 



if^Bf > (|Mr + m|,J(HVml,J, 
2B^i < 2|^|Vm2 +mjj. 



(7) 



In fact one should minimize the full 1-loop Higgs po- 
tential, but as the first approximation we consider only 
the tree-level. We compensate this by adding radia- 
tive corrections, which contain functions of particle mass 
eigenstates generated by elektroweak symmetry breaking 
(EWSB) mixing. In that way a proper EWSB mechanism 
is included in our procedure. 

To obtain the physical masses of supersymmetric par- 
ticles one has to diagonalize proper mass matrices. For 
the down squarks we have 
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m 



+ '^1" " IC^^w +r7i|)cos2/3 



-m<jfc((A_D)fcfe + ^tan/3) \ 
k- + "^di- + U'^w ~ "^D cos 2/3 1 ' 



(8) 



and for sleptons: 



TO^^ - 5(2m^ -m|) cos 2/3 -mefc((A£;)fcfc +^tan/3) 



-mefc((A£;)fefc +^tan^) 



+ + {m^ - m%) cos 2/3 



(9) 



where = cos^ 9w , and the Weinberg weak mix- 
ing angle is sin^ 9w — 0.22. We have denoted by m^ik 
and TOgfc, k = 1,2,3, the masses of down quarks and 
charged leptons respectively. The L-handed elements 
are the eigenvalues of the running mass parameters iiIq 
and whereas the i?-handed ones are the eigenval- 
ues of the singlet parameters m|, and m|;. They are 
obtained from the RGE procedure after diagonalization. 
The diagonalization procedure involves multiplication by 
orthogonal matrices which in the standard trigonomet- 
ric parametrization introduce mixing angles between the 
weak and mass eigenstates in the following way: 

= di cos + sin6', 
dj^ = —dis'md + d2 cos9, (10) 

cl = ei COS0 -f 62 sin0, 

Sr = — ei sin -t- 62 cos 0. (11) 

For completeness we list below the explicit expressions 
for the mixing angles and mass eigenstates. For the d- 
type squarks we finish with 

sin26''' = -2m^fc((A£i)fefc /itan/?) 

X [(m|, - m|, - 0.34m| cos 2(3)^ 

+ 4m2,((Az))fcfc-)-/itan/?)2]-i/2^ (12) 

{Pi.D)kk being the fc-th eigenvalue of A^j, and 

'"dl, = ^Kl+'^J5^)-^"4cos2/3 

/ (A£,)fcfc +pitan/3\ 
+ raa.[m,.T — )■ (13) 

For sleptons the analogous expressions read: 

sin20'^ — — 2mgfc((A£;)fefe + /itan/3) 

X [(mgfc -mlk -0.04to| cos 2/3)2 

+ Aml^{{AE)lk+^itanPf]-^'^ (14) 

2 1/2, 2\ ^2 n r, 

^^1,2 2 4 

/ (A£;)fefc -f ^tan/3\ 

+ mf,k\mekT ^^r-T ■ (15) 

V sm 26'^ / 



Each obtained solution is checked against various con- 
ditions. These are (1) finite values of Yukawa couplings at 
the GUT scale; (2) requirement of physically acceptable 
mass eigenvalues at low energies; (3) FCNC phenomenol- 
ogy (6 — > S7 processes). 

The model we have described has only five free param- 
eters: three GUT parameters Aq, toq: ^i/2, the sign of 
/Lt and tan /3. 

III. NEUTRINO MASSES AND TRANSITION 
MAGNETIC MOMENTS 

Many different proposals of generating neutrino mass 
matrix may be found in the literature 19]. In the present 
paper we use the loop mechanism of generating Majorana 
neutrino masses, possible in RpV models |9l llOl IllL ll^ . 
It is also well known that by considering photon interac- 
tions with the particles in the loop a transition magnetic 
moment is generated. 

A. The neutrino masses in R-parity breaking 
MSSM 

Let us first recall the known results. Within the R- 
parity breaking MSSM neutrinos are, in general, massive. 

In the lowest order, the contribution to the neutrino 
mass matrix reads 

Mtr = A,A,, gl 

Ml + M2 tan^ Ow 
^ 4(^m2j,(Mi + M2 tan^ Ow) sin 2/3 - MiM2iJ?) ' 

(16) 

where = — VdKi, and (i^i) are the vacuum expec- 
tation values of the sneutrino fields. 

Going beyond the tree level one may consider diagrams 
depicted on Fig. ^ The down-squark mixing is realized 
through the angle 9. The quark mixing may be taken into 
account in two ways: as mixing of down-type quarks and 
as mixing of up- type quarks (see discussion in Ref. 20j). 
Since in our case the d-quarks enter the loops, only their 
mixing will infiuence the results. The case of w-quarks 
mixing is equivalent to switching the mixing off and has 
no impact on the outcome of our calculations. We denote 
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FIG. 1: Feynman diagrams representing the squark-quark 
loop contribution to the Majorana neutrino mass. 
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FIG. 2: Feynman diagrams representing the slepton-lepton 
loop contribution to the Majorana neutrino mass. 



by V the standard Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix. 

Altogether the contribution to the neutrino mass ma- 
trix coming from the squark-quark loop may be expressed 
as 



jkl I \ a J 



(17) 



where the loop integral is 



sin 261* 



In T 



jk 



In X 



Jk 



"jk 



1 



,,jk 



1 



(18) 



In the above expressions we have introduced two dimen- 



sionless quantities, 



Jk 



-di 



/m|fc and xl^ 



jk 



m 



do 



The factor 3 in Eq. (|17|l comes from the summation over 
three quark colors. 

The situation is simpler for the slepton-lepton loop 
(Fig. 121), as the lepton mixing is much weaker and we 
feel justified to neglect it. By analogy with Eq. ((T7jl we 
obtain 



1 



167r2 E 

jk 



>^i]k>^t'kj{vjk'meJ +vljm^k), (19) 



with the loop integral taking the form 

^ sm2(j)^ I Inj/^*" Iny^'" 

where now vf'-^ = m^^j j rr?-k and = rr^^jjrr?-, 




(20) 



The functions wj^. and 



•j j,, which are free of i?-parity 
breaking SUSY paprameters, can be calculated by run- 
ning the MSSM RGE, in which the whole low energy 
particle/sparticle spectrum is generated. In our calcu- 
lation we have used the following values for the quark 
sector: m„ = 5 MeV, = 9 MeV, = 175 MeV, 
= 1.5 GeV, mb — 5 GeV, mj = 174 GeV, as well as 
the CKM matrix in the form 



0.9755 0.2195 0.003l\ 
V = I -0.2201 0.9746 0.0390 
-0.0055 -0.0387 0.9989/ 



(21) 



This corresponds to the sines of quark mixing angles 
si2 = 0.2195, S23 = 0.039, and S13 = 0.0031 20] in 
the standard trigonometric parametrization. We have 
neglected the possible appearance of the CP violating 
phase in the quark sector. 

The values of the calculated loop integrals have been 
arranged for clarity in the form of matrices. For the in- 
put parameters ^0 = 100 GeV, toq — mi^2 — 150 GeV, 
tan/? = 19, /X > we get: 



V 



2.361 X lO-"* 
« - I 2.361 X 10-4 
2.358 X lO-'' 



4.593 X 10-3 
4.593 X 10-3 
4.587 X 10-3 



V — 



^6.547 X 10 
6.547 X 10 
16.541 X 10 



1.310 X 10 
1.310 X 10 
1.308 X 10 



1.687 X 10-1^ 
1.687 X 10-1 
1.684 X 10-1 

2.528 X 10-1' 
2.528 X 10-1 
2.525 X 10-1 



For the second set of parameters Aq = 500 GeV, mg 
TOi/2 = 1000 GeV, tan/3 = 19, fi > we find: 



^3.396 X 10-5 
3.396 X 10-5 
,3.396 X 10-5 



6.606 X 10-4 
6.606 X 10-4 
6.605 X 10-4 



2.419 X 10-2^ 
2.419 X 10-2 
2.419 X 10-2 
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/1.023 X 10-5 2.046 x lO'^ 3.879 x lO'^X 
/ = 1.023 X 10-5 2.046 x lO'^ 3.879 x lO'^ . 
yi.023 X 10-5 2.046 x 10"^ 3.879 x 10"^ 

B. The phenomenological neutrino mass matrix 

The constraints on the various products of couphng 
constants A^Ai', XijkK'kj, and A^j,A-,j,^- or their val- 
ues can be find by a comparison of theoretical neutrino 
mass matrix calculated within RpV MSSM with the phe- 
nomenological neutrino mass matrix, derived from anal- 
ysis of neutrino data by making viable assumptions. It 
is usually done by assuming that different contributions 
to theoretical neutrino mass matrix do not significantly 
compensate each other, i.e., it is possible to extract limits 
on individual contributions without knowing the others. 

The neutrino phenomenological mass matrix in flavor 



basis, Ai^^, can be written as 

= U* M'^'^a U\ (22) 

where 

j^^diag ^ diag{mi, ma, mg). (23) 

Here are the moduli of neutrino mass eigenvalues. 
The three neutrino mixing scheme is assumed, namely 

i^aL = VaiViL-, a ^ e,^,T, i ^ 1,2,3. (24) 

are flavor neutrino states and Vi are the mass eigen- 
states. 

For massive Majorana neutrinos the unitary 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix 
U in the standard parameterization has the form 



C12C13 S12C13 sise 



U = -S12C23 - ci2S23Si3e** C12C23 - si2S23Si3e"^ S23C13 e'"i2 , (25) 



S12S23 - Cl2C23Sl3e** -C12S23 - Si2C23Sl3e"^ C23C13 



1 

























where Sij = sin 9ij , Cij = cos Oij . Three mixing angles 
Oij {i < j) vary between and 7r/2. The 5 is the CP 
violating Dirac phase and ai2, 0113 are CP violating Ma- 
jorana phases. Their values vary between and 2n. The 
expression for A4^p in terms of rrii, Oij, S, an, a\j, is 
given in the Appendix. 

The analysis of the Super-Kamiokande atmospheric 
neutrino data and the global analysis of the data of 
the solar neutrino experiments and KamLAND experi- 
ment ^ yield the following best fit values of the relevant 
neutrino oscillation parameters: 

Am^3 = 2.1x10"^ eV^ sin^ 26I23 = 1.00, 
b.m\^ = 7.1x10^5 gv2^ tan^^ia = 0.40. 

(26) 

The neutrino mass-squared difference is determined as 
Amfj, = m| — mf. For the angle 6*13 only the upper 
bound is known. From exclusion plot obtained from the 
data of the reactor experiment CHOOZ '2l'| we have 

sin2 6'i3 < 5 X 10^2 (90% c.L). (27) 

The CP violating phases (5, a\i and 0113 remain undeter- 
mined. 

Neutrino oscillation data are insensitive to the absolute 
scale of neutrino masses. The values of neutrino masses 
depend on the lightest neutrino mass, on the neutrino 
mass spectrum and the neutrino mass-squared differences 
AmJa and Amjg . The neutrino oscillation data are com- 
patible with two types of neutrino mass spectra: 



i) The normal hierarchy (NH) mass spectrum, which cor- 
responds to the case 

TOi ^ m2 < 7713. (28) 

ii) The inverted hierarchy (IH) of neutrino masses. It is 
given by the condition 

777.3 ^ < (29) 

The Oj^/3/?-decay is one of the most sensitive known 
ways to probe the absolute values of neutrino masses 
and the type of the spectrum. The most stringent 
lower bound on the half-life of 0://3/?-decay were ob- 
tained in the Heidelberg-Moscow ^^Ge experiment p3 | 
|.yCk.-ea:p ^ Y. yr). By assuming the nuclear 
matrix element of Ref. we end up with |to^/3| = 
Vl^ mi + U^2 "^2 + U^3 W3 < 0.55 eV, where U is the 
neutrino mixing matrix Eq. H25|l . With this input limit 
we can find the maximal allowed values for the matrix 
elements Ai^^ of the neutrino mass matrix, which are as 
follows: 

/O.55 1.29 1.29 \ 
^j^ph-HM^ < 1.29 1.35 1.04 eV. (30) 
\ 1.29 1.04 1.35 J 

In the calculation we used best fit values of mass-squared 
differences ATnfa, A7TI23 and mixing angles 6*12, 6*23 0, 
0: sin2 6'i3 = [0, 5 x lO^^], S, ai2,ai3 = [0, 27r] and 
the whole allowed mass parameter space of neutrinos was 
taken into account. 
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The phenomenological neutrino mass matrix can be 
constructed from the assumption of normal or inverted 
hierarchy of neutrino masses. If mass squared-differences 
and mixing angles of H2t)l) and (|27|) are considered we 
obtain 



(1.80,4.52) 
(0.025,3.12) 
(0.025,3.12) 



eV 

(0.025,3.12) (0.025,3.12) 

(0.028,2.40) (0.945,2.39) 

(0.945,2.39) (0.028,2.40) 

eV 



^j^ph^NH^ = 10- 

/ (0.53, 45.2) (26.9,98.9) (26.9,98.9) 
X (26.9,98.9) (173,254) (182,254) 
\ (26.9,98.9) (182,254) (173,254) 

Here, we assumed that the mass of the lightest neutrino 
is negligible. The allowed ranges have origin in the un- 
certainty coming from the ^13 parameter and the CP vi- 
olating phases. 



C. The magnetic moments of neutrinos 

Once neutrinos are massive, they can have magnetic 
moments. We distinguish magnetic moments of Dirac 
and Majorana neutrinos: 

i) The Dirac magnetic moment^ which connects left- 
handed electroweak doublet neutrino to a right- 
handed electroweak singlet neutrino Vjji {i, j — e, /i, r). 
One can express the effective Hamiltonian, which con- 
serves the total lepton number, as 

HSf - l^^u^,>^^L<T''^'^JRFpa + h.c, (31) 

where /i^o is Dirac diagonal {i ~ j) or transition {i ^ j) 
magnetic moment between states and I'jn. 

A minimal extension of the SM with a right-handed 
singlet neutrino yields a diagonal neutrino magnetic mo- 
ment of El 



3 Gprnem^ 



47r2 



V2 



Here, fJ-B — e/(2me) is the Bohr magneton and is the 
neutrino mass. By using the neutrino oscillation data 
one finds > 4 x 10"^° hb- It is believed that larger 
values of magnetic moment are possible in extensions of 
the SM, e.g., in models with large extra dimensions [25l |. 
A neutrino magnetic moment /i^ of the order of 10~^^ /is 
might be an explanation of solar neutri no p roblem by 
flipping j/e to sterile neutrinos [2II27L I2II29I I3(ll3l|| . 

ii) The Majorana magnetic moment acts between neu- 
trino fields of the same chirality, namely i^iL and v^j^ As 
a consequence it violates the total lepton number by two 
units (AL = 2). The effective Hamiltonian takes the 
form 



1 



a/3, 



h.C. 



(33) 



Majorana neutrinos cannot possess a flavour diagonal 
{i = j) magnetic moment due to the CPT theorem. They 
can have only transition {i ^ j) magnetic moment j32j . 
From the deflnition one sees that fi^. . vanishes for i — j 
and also that magnetic moment is antisymmetric in in- 
dices i,j: IJLy^. = ~lJ.yj,. 

The limits on neutrino transition magnetic moments 
arise from laboratory measurements of the Ve — c scat- 
tering cross section using solar, atmospheric and terres- 
trial neutrinos |33| . These experiments have put upper 
bounds on jj,^, . , which are as follows: 



A^f^^ < 6.8 X 10-1" fiB, 
fi^^ < 3.9 X 10-^ t,B. 



(34) 



The most restrictive limits remains those from the atro- 
physical analysis of the energy-loss rate of stars due to 
the 7 — > z^77 process, which are j34j 



mS, < 3 X 10 



-12 



Mb- 



(35) 



The assumption is that this energy loss mechanism of 
giant red stars cannot exceed energy loss via weak pro- 
cesses. For Majorana transition magnetic moments /i*^^ 
the above limits are twice stronger j35l |. 

The cases of Dirac and Majorana magnetic moments 
have fundamentally different physical applications and 
need to be considered separately. In this paper we as- 
sume massive neutrinos to be Majorana particles. In ad- 
dition, we shall take advantage of the fact that neutrino 
magnetic moments have an intimate connection to their 
masses, which is model dependent. Recently, the contri- 
bution of transition magnetic moments to the neutrino 
mass matrix was estimated in ,,35j . It was found that if 
neutrino transition magnetic moments are of the order 
of their current upper bound, their contribution to Mij 

(32) 

can exceed the experimental value of A^fj*- The pur- 
pose of this paper is to discuss the impact of the -Mfj 
on the transition magnetic moments of Majorana neu- 
trinos within the GUT constrained MSSM with explicit 
i?-parity violation. 

In order to generate the magnetic moments within the 
i?-parity breaking SUSY one needs to attach a photon 
to one of the internal lines in the loops from Figs. Q and 
121 For each Feynman diagram we have two possibilities. 
One can attach the photon to quark (lepton) or squark 
(slepton) internal line of the loop. However, the con- 
tributions coming from photon-sparticle interactions are 
suppressed by big masses of the latter and we may safely 
ignore them. 

The contribution coming from the squark-quark loop 
takes again into account squark and quark mixing. We 
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end up with 



we obtain 



(36) 



where now the loop integral takes a slightly more com- 
plicated form 



q 



sin 261'' 



jk T jk 



1 



(1 



(X2 ^ Xl) 



(37) 

Here Qd — 1/3 is the d-quark charge in units of e, and 
TOg denotes the electron mass. We work in mathematical 
units, ie. we set h = 1. We note that the resulting 
formula is antisymmetric as expected. 

The contribution from the slepton-lepton loop reads 



4QeTOe 

167r2 



(38) 



/ , '^ijk'^i'kj 



where the loop integral is as previously equal to 



w 



sin 2(j)^ 



jk 



' jk 1 jk 



jk 
V2 



y'2? 



{V2 



yi) 



It is worthwhile to notice that if quark mixing is 
neglected, sin2(?''^/2 « m^^fc m^fc /(m^j. — "^|fc) and the 

masses of squarks m^k , m^k are replaced with their av- 
erage value TOjfc, the expression for fil ^ of Ref. is 
reproduced. The same procedure allows to obtain the 
result of 23 also for fJ,^^,,,- In this case sin20''/2 « 
m^knigk / {m^k — rn^k) and rrigk, m^k = m^k (the aver- 
age mass of sleptons) have to be considered. 

The matrices w^j. and wjj. in flavor space are evaluated 
for two representative sets of input parapeters Aq, mo, 
mi/2, tan/3 and /i > 0. By assuming Aq = 100 GeV, 
mo = mi/2 = 150 GeV, tan/3 = 19, /i > we end up 
with 

/2.582 X 10-12 5.024 x 10-" 2.332 x 10"' 



6.759 X 10-1° 
y 2.750 X 10-^ 



1.315 X 10-^ 6.068 X 10~ 
5.350 X 10-^ 2.421 x lO"^ 



^1.223 x 10-1* 2.446 x 10-12 5.163x10-"' 
= I 2.691 X 10-1" 5.384 x 10-** 1.132 x 10-^ 
^4.832 x 10-^ 9.668 x 10-6 2.019x10-* 

If larger SUSY mass scale is assumed with Aq = 
500 GeV, mo = mi/2 = 1000 GeV, tan/3 = 19, /x > 0, 



/9.895 X 10-15 1.925 x 10-13 8.809 x lO-i^^ 



2.780 X 10-12 
y 1.383 X 10-^ 



5.408 X 10-11 
2.690 X 10-s 



2.466 X 10-3 
1.216 X 10-6 



w 



( 5.316 X 10-1^ 1.064 X lO-i* 2.126 x lO-i^^ 

1.301 X 10-12 2.602 X 10-1" 5.195 x 10-^ 
y2.755 X 10-1" 5.511 x 10"^ 1.098 x 10-^ 



IV. RESULTS 

The main purpose of our present work is to calculate 
the transition magnetic moments for Majorana neutri- 
nos. To achieve this goal one proceeds as follows. First 
one finds the values of the loop integrals w^, w*, and 
within some GUT scenario. We have chosen two sets 
of parameters: (Aq = 100 GeV, toq = "1-1/2 = 150 GeV) 
and {Aq = 500 GeV,TOo = ™i/2 = 1000 GeV), in both 
of them keeping tan /3 = 19 and positive /i. This allows 
us to construct the theoretical mass matrices, Eqs. (|17|) 
and 1)19(1 . and confront them with phcnomcnological ones. 
Next we calculate the contributions from the ^-loop and 
the (?-loop. The crucial point is that we consider each 
mechanism separately, which means that only one ele- 
ment from the sums in Eqs. H17|l and H19|l is picked up at 
a time. Such approach is justified by the usual assump- 
tion that there is no fine-tuning between different con- 
tributions, which therefore may by analyzed separately. 
Explicitly, the simplified expressions read: 



(l-<5..')A^?.i^USY> 



^ViaVlawUmq. 

a 



(40) 



^ (l-<5.,0-M?.'/susY, 



(41) 



where the [...]max symbol denotes that we have checked 
all the combinations of indices {j, k, V\ in the first case 
and {j, fc} in the second one, and picked up the dominant 
mechanism. 

The coefficients /sugy convert the neutrino masses 
into magnetic moments. We list their values, obtained 
for four different sets of GUT parameters, in Tab. HJ 
The /gusY^*^ parameter takes into account quark mix- 
ing, which is absent in /gugy This additional mechanism 
slightly lowers the maximal allowed value of the magnetic 
moment One sees also that the dependence on the tan /3 
parameter is very weak, therefore in the following we will 
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TABLE L The SUSY coefficients converting eiements of neutrino mass mass matrix to transition magnetic moments of 
Majorana neutrinos (see Eqs. 1401 and 14111 for definition), ft is assumed > 0. 





SUSY input 








1 fie SUSY conversion 


coefncient 








mo 




tan P 


f 

JSUSY 




rq-CKM 
/sUSY 




/sUSY 




[GeV] 


[GeV] 


[GeV] 




[eV-1] 




[ eV-^] 




[ eV-1] 




100 


150 


150 


5 


(0.3, 1.0) X 10" 


-16 


(2.8,8.8) X 10" 


-17 


(0.5,1.5) X 10" 


-15 








19 


(0.3,1.2) X 10" 


-16 


(2.8,9.8) X 10" 


-17 


(0.5,1.6) X 10" 


■15 


500 


1000 


1000 


5 


(1.1,2.8) X 10^ 


-18 


(1.0,2.4) X 10" 


-18 


(1.7,4.1) X 10" 


-17 








19 


(1.1,3.1) X 10" 


-18 


(1.0,2.7) X 10" 


-18 


(1.7,4.3) X 10" 


-17 



TABLE If: Ttie Majorana neutrino transition magnetic moment (ij — efi, er, fir) derived from phenomenoiogicai mass 
matrices (see text for detaiis). Ttie upper bounds on magnetic moments are obtained by using ttie current neutrino osciiiation 
and tfie Oz^/3/3-decay data. Tfie predictions for fi^.^ were caicuiated by assuming the inverted hierarchy or the normai hierarchy 
of neutrino masses and that the iightest neutrino is massiess. The ranges of aliowed vaiues correspond to uncertainty in SUSY 
coefficient fg'ljgY ^'^'^ CP-vioiating vioiating phases. Two different GUT scenarios are considered. The remaining SUGRA 
parameters are tan/3 = 19 and /i > 0. /is is the Bohr magneton. 





SUSY input in GeV 




Transition mag 


;netic moment /i^.^ in ^lb 






mo 


mi/2 


ij 


QvPP constraints 


inverted hierarchy 


normai hierarchy 










iept' 


on-slepton ioop mechanism 




100 


150 


150 


ejx, er 


< (0.62,2.10) X 10"^^ 


(0.01,5.00) X 10"" 


(0.13,1.60) X 10"" 








HT 


< (0.50, 1.70) X 10"^^ 


(0.45,3.80) X 10"" 


(0.87,4.10) X 10"" 


500 


1000 


1000 


ei-t, er 


< (2.20,5.50) X 10"^^ 


(0.04, 13.0) X lO"^'' 


(0.47,4.20) X 10"^^ 








HT 


< (1.80,4.50) X 10"^'^ 


(0.16, 1.00) X 10"^* 


(0.32,1.10) X 10"^* 










quark-squark loop mechanism (without d-quarks mixing) 


100 


150 


150 


e/i, er 


< (0.41, 1.50) X 10~^® 


(0.08,36.0) X 10"" 


(0.08,1.10) X 10"^* 








/ir 


< (0.33,1.20) X lO"^'^ 


(0.30, 2.80) X 10"^* 


(0.58,3.00) X 10"^* 


500 


1000 


1000 


e/i, er 


< (1.40,4.00) X 10"^* 


(0.03,9.60) X 10"^" 


(0.30,3.00) X 10"2° 










< (1.20,3.20) X 10"^** 


(1.00,7.30) X 10"^" 


(2.00, 7.80) X 10"^° 










quark-squark ioop mechanism (with d-quarks mixing) 


100 


150 


150 


e/i, er 


< (0.36, 1.30) X lO-^'^ 


(0.07,31.0) X 10"^^ 


(0.76,9.70) X 10"^^ 










< (0.29, 1.00) X 10"" 


(0.27,2.30) X 10"^* 


(0.51,2.50) X 10"^* 


500 


1000 


1000 


e/i, er 


< (1.30,3.40) X 10"^** 


(0.03, 8.30) X 10"^" 


(0.28,2.60) X 10"^° 








/ir 


< (1.10,2.80) X 10"^** 


(0.97, 6.40) X 10"^'' 


(1.90,6.80) X 10"^° 



discuss only the case of large tan/3 = 19. For complete- 
ness we include also in the table the lower bounds on the 
/ coefficients, which were obtained by finding the small- 
est of the expressions in square brackets in Eqs. (|40|l and 

The resulting values of /i^ ., for six different scenar- 
ios are presented in Tab. ^ Due to antisymmetricity 
the diagonal elements are zero. The ranges of val- 
ues in the first column are the highest upper limit and 
lowest upper limit on the transition magnetic moment, 
coming from different mechanisms (different combination 
of coupling constants AA or A'A'). They were calculated 
using the J^p^^^^^ mass matrix. One clearly sees that 
in this case the different mechanisms give comparable 
results. However, the imposed GUT constraints have a 
much bigger impact and introduce two orders of mag- 
nitude differences. The values from the first column 



for the first set of GUT parameters {Aq — 100 GeV, 
Too — mi/2 — 150 GeV) may be roughly compared with 
the predictions published in Ref. [l^ showing that for 
this special case we our results are compatible with the 
previously published. 

The two remaining columns present the ranges of al- 
lowed values of the magnetic moments if one assumes nor- 
mal or inverted mass hierarchy of the neutrinos. These 
results take into account various mechanisms of generat- 
ing /ii^. ., as well as the experimental uncertainties shown 
in mass matrices M^'^~^^ and Mp'^^^^ , In this case 
the ranges span roughly over 1-2 orders of magnitude. 
It is visible that, regardless of the GUT parameters, one 
should not expect the magnetic moment to be greater 
than 10~^''/iB, which is not possible to detect experi- 
mentally nowadays. 

It is interesting to investigate the infiuence of quark 
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10 ^ 



1= I I I mill — I I 1 1 iiiij — I I 1 1 iiiij — I I Mi 



IH, ij=e|i, ex 
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10 = 



10-^^ k 



lO^g ^3 [eV] 

10 N I I I llllll 1 I I I Mill 1 I I I Mill 1 I I l> 



NH, ij=e|i, ex 



10-^^ k 




= 111 lllllj 1 1 1 Mill 

= IH, ij=|ix 


1 1 1 Illlj 1 1 LWItt 
— 




1 


1 



m3 [eV] 



till lllllj — I I I llllll — I I 1 1 Mill — I I \y\m 



XQ" 10^ 10 
[eV] 



lo^io-' 




I I 



10"^ 10^ 10 
[eV] 



FIG. 3: Transition magnetic moments of Majorana neutrinos /i^^^ , ^i,^^ and ^.v^^ as function of the lightest neutrino mass ma 
(the inverted hierarchy of neutrino masses, upper panels) and mi (the normal hierarchy of neutrino masses, lower panels). The 
results are presented for sin^ ^13 — (the region with solid line boundaries) and sin'^ 613 = 5 x 10~^ (the region with dashed 
line boundaries). In the calculation the best fit values of mass-squared differences Amf2, Am23 and mixing angles 612, O23 were 
considered For CP- violating phase we assumed S, ai2,ai3 = [0,27r]. The largest value of the SUSY coefficient fsusY 

for the set of SUSY parameters {Ao — 100 GeV, mo = mi/2 = 150 GeV), tan/3 = 19 and > was taken into account. IH 
and NH denote the inverted hierarchy and the normal hierarchy of neutrino masses, respectively, is the Bohr magneton. 



mixing on the results. When quark mixing is switched off, 
the maximal allowed fj^l values drop down approximately 
by 10 to 20%. 

In the case when the lightest neutrino mass is not set to 
zero, one may express all the remaining mass eigenstates 
in terms of rriK^^) and appropiate differences of masses 
squared, which are known from the results of the neutrino 
oscillation experiments. On Fig.|31we present the allowed 
values for different elements of the neutrino mass matrix 
as the function of mi(3) in the case of NH and IH. The 
shaded regions correspond to all possible combinations of 
the two Majorana ai2 and ai3, and one Dirac phase S, for 
sin^ 013 = (dark shading) and sin^ 9i3 < 0.05 (Hghter 
shading). One sees clearly, that the phase factors as well 
as possible non-zero mass of lightest neutrino modify the 
results in a non-trivial fashion. 



V. SUMMARY 

Since the neutrinos have no electric charge they can be 
Majorana fermions. A Majorana neutrino is more nat- 
ural than a Dirac neutrino as it allows to explain the 
smallness of their masses in theories with grand unifica- 
tion. Majorana neutrinos have zero diagonal magnetic 
moments while they may in general have transition mag- 
netic moments fx^,.- {i ^ j). The question is, how large 
they are? In this paper we studied this issue within the 
GUT constrained MSSM with explicit i?-parity violation. 

The transition magnetic moments were induced by at- 
taching photons to the internal lines of lepton-slepton 
and quark-squarks loop diagrams that generate neutrino 
masses, /i^. . has been written as a product of SUSY 
conversion coefficients fsuSY, which reflects the depen- 
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dence on SUSY parameter space and quark mixing, and 
the element of Majorana neutrino mass matrix in flavor 
basis. 

Two possible GUT scenarios were considered: = 
100 GeV, mo = mi/2 = 150 GeV, and = 500 GeV, 
Too ~ TOi/2 ~ 1000 GcV, in both of them having 
tan /? = 19 and positive /i parameter. Our study showed 
that fsuSY is nearly insensitive on the value of tan/3 
and that its dependence on the scale of GUT common 
mass and coupling constants parameters mo, mi/2, and 
Aq dominates. The difference in fsuSY for these two 
scenarios reach around two orders of magnitude. It is 
worth mentioning that contributions to /i^. from lepton- 
slepton mechanisms dominates over contributions from 
quark-squark diagrams. 

The transition magnetic moments were calculated by 
considering three differently constructed phenomenolog- 
ical neutrino mass matrices. The elements of the first 
one (yV/JP'*"^*'^) represent the maximal allowed values, if 
neutrino oscillation data 0, 0| and the lower bound on 
the Qvfijd-Aec&y half- life of ''^Ge jl^l are taken into ac- 
count. The corresponding upper limits on ^^^^ , /i^^^ and 
/ijy^^ have been found to be of the order of 10"^*^ hb and 
10^^^ /is for the assumed two SUSY scenarios. 

Two phenomenological neutrino mass matrices were 
generated by considering normal and inverted hierarchy 
of neutrino masses and that the mass of the lightest neu- 
trino is zero. This allowed to make predictions for tran- 
sition magnetic moments of Majorana neutrinos. For 
low (large) SUSY parameter scale we found the maximal 
value 5 x 10^^'' (1-3 x 10^^^ /i_b). The dependence of 
the results on the neutrino mass pattern is weak. 

The magnetic moments /i^^. ^ , and ^ were cal- 
culated also as a function of the lightest neutrino mass. 



The best fit values of neutrino oscillation parameters were 
considered. Only for the mixing angle ^13 it was assumed 
sin^ 6*13 = and sin^ 6'i3 = 5 x 10"^. It is worth mention- 
ing that the maximal mixing of atmospheric neutrinos 
(^'23 — T^/^) implies that the allowed ranges of ^j^^^ and 
/ijy^^ are the same, if the Dirac and Majorana CP vio- 
lating phases vary from to 7r/2. The results show that 
especially ^iv^^ ex (inverted hierachy of neutrino masses) 
depends strongly on the values of CP phases. 

The obtained results show that the values of the Majo- 
rana transition magnetic moments might be significantly 
above the scale of the Dirac-type magnetic moment in 
minimal extension of the SM with right-handed neutri- 
nos. However, the maximal values are still too small to be 
tested in the present laboratory experiments or to have 
astrophysical consequences. 
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VI. APPENDIX 

In this appendix, for covenience of the reader we 
present explicit expressions for the neutrino mass ma- 
trix in flavor basis, M.af} {a,/3 = e,/j,,r), as functions 
of the moduli of neutrino mass eigenvalues rrii , of mixing 
angles 0ij and of the CP violating phases S, ai2, ais [36j. 
The matrix element M is symmetric, i.e., it is defined by 
six elements: 
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Mee = Ci3Ci2mi + C^g S^j m2 6 '^"12 + S^g 6^ ma 6 '^"13, 
A^e;* = -C12C13 (C23S12 + C12 S23 Sl3 e""* ) TOi 

+ Ci3 S12 ( C23 C12 - S23 S12 si3 e"**^) m2e"'^"i2 + C13 S23 S13 e**^ m3 e"*^"" , 

Mer = - C12 Ci3 (- S23 S12 + C23 C12 Si3 6"'*) mi 

-C13S12 (ci2S23 + C23 S12 si3 e"'* ) m2e"'^"i2 + C23 C13 S13 e'* ms e"'^"i3 , 

■^W = (C23 S12 + 2 C23 C12 §23 S12 Sl3 e"*'' + C12 S23 Sl3 e"^*^) TOl 

+ (c23Ci2 - 2 C23 C12 S23 S12 si3 e"*'' + sl^sl^sl^e-^^^) mae-^^^is + c^g s^g m3 e-'^^ia , 

A^/HT = - (C23 S23 S?2 - C23 C12 S12 Si3 e~'* + C12 S23 S12 Sl3 6""' - C23 S23 ^13 e~^"^ ) mi 

- (c23 C12 S23 + C23 C12 S12 si3 e~"' - C12 S23 S12 si3 e~'* - C23 S23 S12 Si3 e~^"^ ) m2 e~'^"i2 
+ C23 Ci3 S23 ms e"*^"i3 , 

■A^TT = (.S23SJ2 - 2c23Ci2S23Sl2Sl3e""' + C23Ci2S?3e"^*'^)TOi 

+ (c?2 S23 + 2 C23 C12 S23 si2 S13 c"*'^ + C23 S12 s^g e-^*-^) m2 e"*2"i2 + C23 TO3 e-*2"i3 . 

(42) 



The maximal mixing of atmospheric neutrinos (^23 = 
7r/4) impHes symmetry among some of elements of the 
neutrino mass matrix Aiafs- If CP-violating phases are 
considered as free parameters the allowed ranges of lAlgjul 
and lA^erl (l-^^pl and lA^rrl) coincide each with other. 
If additional restriction is taken into account, namely 
sin2(0i3) = 0, we obtain 

Mee = cf2mi + Si2'^2e~*^"i2 

v2 

Mer = --y=Ci2Si2 (- + e"'^"i2 ) , 

■^w = ^ (s?2 TOl + c?2 m2e-'2"i2 + m3e-'2«i3) , 

Mi^r = - ^ ( s?2 + cf2 m2 e-'2"i2 - m3 e-'2"i3 ) , 

Mrr = ^ (512 TOl + C12 m2e~*^"i2 + m3e~'^"i3) . 

(43) 



In this case IA^bmI = l-^erl and A^/a/i = A^rr- In addi- 
tion, for 5, ai2, ai3 between and 27r the maximal and 
minimal absolute values of all three elements Ad^r, -M-i^in 
and Mtt are the same. 
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